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ABSTRACT 



An investigation into the problems of economically pro- 
cessing sewage on board Naval ships resulted in the develop- 
ment of two computer simulations employing Monte Carlo 
analysis to describe the generation of sewage. Simulation 
one was based on a non-homogeneous Poisson process. For 
simulation two, an empirical distribution describing the 
arrival behavior of sewage to the holding/processing unit 
over a 24 hour period was applied to known data on sewage 
generation . 

Results of the two simulations were compatible with one 
another. Aside from pointing out a most feasible combina- 
tion of holding tank capacity, processor rate and processing 
policy, the simulations also indicated that a revision of 
the Navy's design parameter for the daily per capita sewage 
generation rate was in order. The simulations were designed 
so that use by others with their own data would be an easy 
matter . 
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I. 



INTRODUCTION 



A . BACKGROUND 

During the time between the conception of this thesis 
and its completion, significant changes in both policy and 
attitude have been made by the Navy in the field of pollu- 
tion abatement. These alterations have resulted from in- 
creasing public pressure for the- ecological restoration of 
U.S. inland waters and subsequent congressional legislation 
which outlined future restrictions to all seagoing vessels 
and watercraft. More specifically, in August 1971, the 
Chief of Naval Operations stated "that the water pollution 
abatement program for Naval ships has as its long term goal 
the installation of two basic systems aboard ship: 

1. A Marine Sanitation Device (MSD) which processes 
sewage to an acceptable effluent; 

2. A Collection, Holding and Transfer (CHT) System, 
which is designed to transfer all liquid wastes (except 
oil) ashore or to a barge." 

Late in 1971, however, the Senate issued further re- 
strictions to vessels to the extent that, by the year 1981, 
there shall be no toxic effluent discharge, whatsoever, 
from any ship into U.S. inland waters. As a result. Naval 
design efforts refocused primarily on the development of 
a Collection, Holding and Transfer System. This shift in 
emphasis was also aided by the fact that prototypes of 
Marine Sanitation Devices had suffered from low system 
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reliability and excessive costs. The present goal, then, 
of the Navy is the development of a sewage holding tank 
system with the capacity to retain 12 hours of sewage gen- 
erated by a full complement crew and with the capability of 
transferring such sewage to an external processing facility 
ashore. Such a system is to be installed in all ships re- 
maining in the fleet after 1977. 

B. STATEMENT OF THE PROBLEM 

It is the contention of the author that neither strictly 
holding sewage for later transfer, nor immediately process- 
ing the sewage internally as it arrives to the processing 
unit, is the most feasible economic solution to pollution 
abatement. A policy of holding sewage and transferring it 
to some other processing station would greatly hamper the 
Navy's capability to act as a mobile self-sustained unit; 
likewise, the policy of rapid processing aboard would re- 
quire expensive units of high output and reliability. It 
is believed that viable trade-offs can be realized through 
a combination of the two, employing a smaller, cheaper pro- 
cessor which could operate at a lower processing rate, and 
a smaller holding tank that could absorb sewage flow during 
peak periods of the day. 

A second contention, that is perhaps more germaine to 
the current Navy policy of a twelve hour holding system with 
subsequent transfer, is the belief that the presently ac- 
cepted daily per captia rate of 30 gal/man/day of sewage 
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generation is over estimated [Ref. 5]. The pareimeter was 
determined from data gathered by the Marine Engineering 
Laboratory on the daily sewage generation rates per man on 
board four ships. By reasoning that "any treatment unit 
used aboard a Naval vessel must be designed to operate with 
a full complement 24 hours per day," the MEL Report [Ref. 

18] recommended that the rates measured be doubled for de- 
sign purposes. This resulted in a figure of 30 gal/man/day. 
Since the original figures were obtained on a per capita 
basis, and there was little variation, it is felt by the 
author that more accurate results could be obtained by 
applying the unmodified data to a model which also considers 
the traffic pattern of personnel using the facilities over 
a period of time . 

It was the purpose of this study, then, to investigate 
the problem of sewage generation and processing aboard a 
Naval ship, and to develop a model which would best de- 
scribe such a process. This model would be useful in aiding 
decisions on the selection of the most feasible treatment 
system. 

The approach taken, as presented in this thesis, was to 
first formulate the problem analytically and to compare it 
with other research performed in related cases. A computer 
simulation was then developed, motivated by a similar study 
done for the air lines regarding sewage handling on board 
large jumbo jets [Ref. 9]. The simulation, its methodoligies 
and how it can be, used are described in Chapter III. The 
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simulation provides a sensitive model which illustrates the 
trade-offs obtainable through any treatment policy tested. 

It also reinforces the contention that a revision of sewage 
generation rates, combined with closer analysis, can lead 
to savings in holding tank size. 

It should be noted that a simulation is no stronger than 
the basic assumptions made regarding its underlying struc- 
ture, and that such assiamptions should be examined criti- 
cally. With this in mind, the simulations were developed 
using the best information available; however, they were 
also made flexible enough so that, should a designer have 
better arrival data or other policies to test, such informa- 
tion could be easily read into the program to give useful, 
real world results . 
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II . PROBLEM FORMULATION 



A. THE BASIC MODEL 

The overall process of generating sewage over a specific 
length of time, storing it and/or processing it can be de- 
scribed analytically. Further, it can be related to a 
group of similar problems that are generally classified as 
storage problems, examples of which are the collection and 
storage of water in a reservoir subject to some release 
rule, and inventory problems in which demand is random. 

Any such process in which either the input or output is 
random comes under the category of a stochastic process. 

Such is the case with the sewage storage and processing 
problem, in which the input, X(t) , (see Figure 1) into a 
storage tank is random in nature. Here X(t) represents the 
amount of sewage generated by the ship's company at discrete 
points in time, t, corresponding to, say, hours or minutes. 
Two important characteristics of this input are 1) its time 




Figure 1. Schematic of System. 
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dependent arrival behavior to the tank, and 2) its overall 
distribution. It is assumed that X(t) is independent and 
identically distributed for different values of t. For this 
discussion, X(t) is assumed to be solely bodily waste; how- 
ever the model could also be expanded to include "hotel" 
and galley waste. 

- •• -rZ (t) represents the .total .amount of sewage in the stor- 
age tank also at time t. Aside from being a function of the 
input random variable, X(t) , it is also a function of the 
rate, a, at which the sewage is processed, and the particu- 
lar "release rule" or processing policy, ip , that is applied. 
The processing policy is essentially a switching function 
which turns the processor on or off based on either a time 
plan (e.g. on at certain hours of the day) , a tank level 
plan (e.g. turn on when the level of sewage reaches a cer- 
tain height) , or a combination of the two. 

At certain times of the day, the processing rate may 
exceed th^ input, X(t) , while at other times it may lag it. 
If the processing rate is much lower than the input rate, 
Z(t) will behave like a sum of independent random variables. 
For such a case, the tank would be expected to overflow. 

If the processing rate is much higher than the input, then 
the tank may never fill at all, and Z(t) will behave more 
deterministically (or equal zero) . One of the questions to 
be answered, then, is what is the optimal processing rate? 
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In general, Z (t) may be described as follows: 



Z (t) 



K 



if Z (t-l)+X(t) -at>K 



Z( 



t-1) +X ( t) -at f 



if 0<Z (t-1) +X (t) -at<K 



Z (t-l)+X (t) , if a=0 

^0 if Z ( t-1) +X ( t) -at<0 

where K is the capacity of the tank. 

This formulation includes two special cases . The 
is the event that overflow occurs in which case : 



( 1 ) 



first 



Z(t) 



' Z (t-l)+X(t) -at, if Z (t-l)+X(t) -at<K 



( 2 ) 



.K if Z(t-l)+X(t)-at>K 

Here the overflow is Z (t-1) +X (t) -at-K . The second case is 
the event that processing empties the tank completely, so 
that : 



Z(t) 



' Z (t-l)+X(t) -at, if Z (t-1) +X(t) -at>0 

■< 

, 0 if Z ( t-1) +X (t) -at<0 



(3) 



B. OTHER RELATED STUDIES 

Equation (1) through (3) are not unique to the sewage 
generation and collection problem. Rather, they parallel 
those proposed by Moran, Gani , and others [Refs. 15, 8, 10] 
under the broader discipline of the stochastic theory of 
reservoirs or storage, an approach which has attracted con- 
siderable theoretical interest, but has seen less application. 

Initially the theory, as proposed by' Moran, required the 
inflow, X(t), to have a known continuous probability dis*' 
tribution and to be mutually independent of the inflow for 
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any other value of t (all of which have the same distribu- 
tion) . The level in the tank, Z(t) , is subject to some re- 
lease rule, Ip, and the sequence {Z(t+1), t>0} is a Markov 
process, with state transition probabilities that can be 
computed from the inflow probability distribution. Consid- 
erable effort has been expended in finding solutions to the 
steady state distribution of Z(t) for both continuous and 
discrete cases . 

These initial studies were extended by a number of 
authors. One such work by F. M. Lloyd [Ref. 12] showed that 
the theory need not be limited to inflows that were mutually 
independent, but that it could also be applied to serially 
correlated inputs with the results forming a bivariate 
Markov Chain. Ali Khan and Gani [Ref. 1] carried this study 
one step further to consider the time dependent solution 
to the same model. 

An excellent example of the application of the above 
theory which is closely related to the sewage problem is 
the work by Avi-Itzhak and Ben-Tuvia [Ref. 3] . Here the 
authors studied the problem of "how big a reservoir to build 
and what should the pumping capacity be, such that the cost 
per cubic meter of water is a minimum." In order to obtain 
their solution, however, the authors first had to assume an 
approximate distribution for the daily input to the reser- 
voir as well as know the costs involved in operating the 
system. 
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In all of the above cases the underlying requirement is 
that the distribution of X(t) is known or can be approxi- 
mated by a sequence of discrete distributions. 
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III. MODELING THROUGH SIMULATION 



A. JUSTIFICATION OF SIMULATION 

The storage theory approach described above leads to 
perhaps the most precise solution to the sewage storage 
problem. It has, however, certain drawbacks which preclude 
its immediate implementation. The first of these is the 
lack of real world information and data with which to 
derive the necessary matrix of transition probabilities. 
Secondly, the mathematics involved are difficult and cxam- 
bersome . This is especially so when one is attempting to 
solve the problem for different processing policies ("re- 
lease rules") and processing rates, each of which would 

f 

require its own transition matrix. Thirdly, the complexity 
of the distribution of X(t) , the input variable, whose sta- 
tionarity over time is questionable , would require parti- 
tioning the process in order to arrive at a solution. Such 
a non-homogeneous process would only further complicate 
the analysis. Lastly, the problem of a reversing boundary 
brought about by abrupt changes in processing rate, accord- 
ing to the processing policy used, would add even another 
dimension to an already difficult problem. 

Taking the above into consideration, it was decided 
that a computer simulation involving Monte Carlo analysis 
was at least a justifiable, if not preferable, alternative. 
Most important to the justification of a simulation was the 
fact that it could be made amenable to any alterations in 
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processing policies and processing rates. It is also adapt- 
able to adjustments in the input distribution if and when 
additional data on the arrival patterns of sewage inputs 
are obtained. These characteristics have been incorporated 
into the simulations developed in this study. 

B. DESCRIPTION OF THE SIMULATIONS DEVELOPED 

‘ Two computer -simulations were’ 'developed using FORTRAN IV 
for runs on the IBM 360/67 computer. The simulations dif- 
fered only in the method of Monte Carlo analysis used to 
generate the sewage input distribution, X(t), which was then 
applied to the basic mathematical model described in Chapter 
II (equations (1) - (3)). For simulation one, X(t) was de- 
termined by assuming a non-homogeneous Poisson process. For 
simulation two, an empirical distribution describing the 
arrival behavior of sewage to the tank over a 24 hour period 
was applied to known data on sewage generation. The tech- 
niques and mathematics involved are described in subsections 
2 and 3 following an overall description of the basic simu- 
lation (subsection 1) . 

1 . Basic Design Common to Both Simulations 

Aside from the methods used to generate observations 
of the sewage input, the two simulations were identical in 
the way the sewage was processed from the holding tank. 

(See general flow chart of the simulations. Figure 2.) 

First it was assumed that the holding tank used in the 
models had sufficient capacity so that overflow was not a 
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Figure 2. General Flow Chart of the Simulations. 
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problem. This was done in order to get some idea of what 
the maximum sewage level (ZMAX) in the tank would be, ac- 
cording to the processing rate and policy used. This 
further simplified the basic model, so that the only equa- 
tions needed for the simulation were : 



Both simulations involved discrete time steps; for 
model one, the unit of time was 5 minute increments; for 
model two, hourly increments. In both cases the total period 
of time simulated was one week, based on the assumption that 
it was desirable to have any sewage build-up during the week 
processed, and the tank emptied by the end of the weekend. 

A new cycle would then correspond to the start of a new work 
week. It should be noted that each simulation day started 
at 0500 real time. At each discrete time increment, a ran- 
dom observation of sewage input, X(t) , was generated and 
the level, Z (t) , was computed using equation (4) for a 
specific processing rate, a. 



separate subroutines could be called independently and ap- 
plied to the model. The following five policies were used: 




(4) 



Anyone of five processing policies, in the form of 



1) Plan 1. Continuous holding with no processing 
whatsoever (a=0) . 



2) Plan 2. Continuous processing at some rate, a 
to be read in by the programmer. 
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3) Plan 3. Process continuously from 0500 until 
1600 each day; after 1600 stop processing when 
the tank level reaches 50 gal; continue pro- 
cessing at 0500 the next day. For weekends 
stop processing anytime the tank level reaches 
50 gal; continue processing when the tank 
level builds up to 500 gal. (This is a com- 
bination of a level policy, Plan 4, and a time 
policy. Plan 5, as described below. 

4) Plan 4. (A level control policy) Stop pro- 
cessing whenever the tank level goes below 
some, prescribed .limit (50 gallons was used in 
the simulations - the value is read in by the 
programmer) . Continue processing when the 
level exceeds an upper limit (1500 gallons 
was used) . 

5) Plan 5. (A time control policy) Process con- 
tinuously for a fixed portion of each day. 

(For the simulations processing stopped at 
hour 16, 2100 real time and continued at the 
beginning of the next simulated day, 0500 real 
time . ) 

Policies one and two illustrate the two extremes 
possible. The policy of continuous holding (Plan 1) was used 
to determine the amount of sewage that would build up over a 
certain period of time. (The special case of a twelve hour 
holding tank was tested using this policy.) The policy of 
continuous processing (Plan 2) was used to get an indication 
of the minimum tank levels that could be expected over any 
range of processing rates desired. 

The remaining three policies were developed as ex- 
amples of probable real world policies . They also served to 
reflect the response of the simulations to changes in policy, 
for which the simulations produced acceptable results. Of 
the three, policy five would be expected to be the least ex- 
pensive since it simply involves turning the processor off 
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once a day. Policy four, on the other hand, would require a 
level sensing device and some sort of automatic start/stop 
mechanism to the processor and would, therefore, entail a 
considerably larger cost. Policy three would only require 
manual operation and its relative cost would fall somewhere 
in between that of policies four and five. The point to be 
considered for. .all three is the fact that the smaller the 
time required for the processor to operate, the lower the 
overall running cost and the time available for system 
maintenance . 

It is emphasized that the simulations are in no way 
limited to just these policies and that any other policy may 
be tested simply by the addition of a subroutine similar to 
those used and a statement by which to call it. Further 
details on the mechanics of using the computer simulations 
are contained in Appendix D. 

2 . Input Distribution by Simulation One 

One of the most important problems to be solved was 
how to generate reasonably accurate observation of the 
arrivals of sewage to the holding tank. The two parameters 
that had to be considered were 1) the arrival rate and 2) 
the quantity of sewage per arrival . Although a study was 
conducted in 1964 [Ref. 19] to determine the quantity of 
sewage generated per man per day on a Naval vessel, no 
definitive statistics were obtained to describe the arrival 
process over the course of a 24 hour period. Owing to the 
rrature of shipboard routine, which can be considered 
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reasonably standard throughout the Navy, it was evident 
that the arrival distribution would be multimodal over the 
period of the day. With such changes in arrival rates it 
was also clear that the process could not be considered 
stationary for a 24 hour period and that some form of non- 
homogeneous process would be required. 

-In a.related study by Gephart and Balachandran 
[Ref. 9] into the reliability considerations of an airborne 
sewage system for large jet liners, data was collected on 
the arrival behavior of passengers to the restroom facili- 
ties. It was concluded by the authors that the arrival 
process could be considered Poisson over "non-overlapping 
and consecutive time segments." It appeared reasonable, 
then, that a non-homogeneous Poisson process could also be 
applied to the sewage generation problem on board a Naval 
ship . 

In its most general form, the process could also 
be considered a compound Poisson process represented by 

I 

N(t) 

X(t) = Z Yi (5) 

i=l 

where {N (t) , t>0} is Poisson process and {Yj^, n=l,2...} is 
a family of independent, identically distributed random 
variables representing the quantity of sewage per arrival. 
If a distribution for the body waste generated per capita 
per arrival were known, this could easily be incorporated 
into the simulation. For the purpose of simplicity. 



21 



however, the quantity of sewage per arrival was taken to be 
a constant and the compound property was disregarded. Based 
on studies [Ref. 19] which determined the quantity of water 
per flush to be 4.5 gallons, a rather conservative estimate 
of 4.9 gallons per arrival was used. 

The final problem, then, was to determine the best 
estimate, of the.. ar rival „behavior at head facilities over a 
24 hour period. It was decided that the day could logically 
be partitioned into three periods : the first being the 

period between 0500 and 1600 when the majority of the crew 
is on board for the working day; the second being the period 
between "liberty call" and "taps" (1600-2200) and the last 
including the time from 2200 until 0500. The intensity 
function X(t) which describes the arrival rates could then 
be written for weekday periods as : 



( X 1 / 0 <s<t j 

tj<s<t2 (6) 

X3 , t^<sst 

and assuming that the majority of the crew is on liberty 
over the weekend, then a reasonable intensity fianction for 
weekends might be 



X(s') = 



X 2 , 0^s<tj ' 

. X 3 , tj '<s<t 

where Xj and X 3 are identical for both sets of equations. 
Thus the mean value function is ; 
t 

I 



(7) 



m 



(t) =J X(s)ds = Xiti+X 2 (tz-tj )+X 3 (t -t2) 



( 8 ) 
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and 



m(t') = ^ X(s')ds = ^2 1 j ' + X 3 ( t-t j ' ) 

o 

and it can be shown that 



(9) 



e-m(t) [m(t)]'^ 

P{N(t)=n} = , n>0 (10) 

nl 

This then describes a non-homogeneous Poisson pro- 
cess. The values for each were determined by comparing 
the mean value function m(t) , multiplied by 4.9 gallons of 
sewage/arrival, to the expected daily and weekly sewage 
build-up. The details of this are contained in Appendix A. 
The final values used for each X^ are presented in Figure 3. 



h 



Arrival/5 min 



Xi=4.5 



0500 




1600 2200 0500 



Time 



Figure 3. Arrival Rates for Simulation One. ' 

The simulation involved generating Poisson variates 
for each 5 minute increments of simulated time, using the 
above parameters during the periods indicated. The number 
of arrivals was then converted to a quantity of sewage X(t) 
according to : 
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X(t) = N(t) X 4.9 gal/arrival 



( 11 ) 



giving the necessary random input for X(t) which is then 
applied to the basic model as described in section 1 above. 

3 . Input Distribution as Described by Simulation Two 

Although the reasoning behind Simulation One was 
considered sound, it was decided that a second simulation 
could be developed through a different approach. Not only 
would it take full advantage of all known data, but it would 
provide an experiment by which the two simulations could be 
weighed against one another. 

In order to get the best approximation to the true 
arrival distribution over a 24 hour period on board ship, a 
group of experienced Naval personnel were surveyed for their 
estimates of the actual arrival pattern (i.e. hydraulic 
loading) of water closets and urinals over a 24 hour period. 
From such a survey, an empirical distribution was developed 
by employing a technique similar to that used to get job 
time estimates for PERT programming (Program Evaluation and 
Review Technique) . The details of this method of estimation 
and how it was used are contained in Appendix B. The re- 
sults are shown in Figure 4a. 

Figure 4b gives a similar hydraulic loading pattern 
derived for a commercial cargo ship with a 40 man crew 
iRef. 13, p. 58]. Although the routine on a commercial 
cargo vessel is considerably different from that on a 
military vessel, it can be noted that "peak loadings are 
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Figure 4a. Estimation of Distribution of Sewage 
Flow Aboard a Naval Ship. 
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Figure 4b. Distribution of Sewage Flow Aboard 
a Cargo Ship. 
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double the average loading and almost five times the minimum 
loading" for both distributions, and that, in general, they 
compare favorably. 

Once this distribution was established, it was 
necessary to convert it to realistic quantities of sewage 
per hour. Again, the hydraulic shipboard data obtained in 
MEL Report 346/64 [Ref. 19] (which is the basis of current 
Navy design parameters of 30 g.m.p.d.) was used to estimate 
a distribution of the daily per capita sewage flow. This 
data, however, was revised in a thesis by D. B. Campbell 
(a fact not incorporated in design specifications listed in 
Ref. 2, p. E-1) . Combining both sets of information it was 
decided that daily per capita sewage flow rates could best 
be represented by a beta distribution with a most likely 
(mode) value of 14.8 gal/man/day, a maximtun of 20.0 g.m.p.d. 
and a minimum of 11.3.’ As shown in Appendix C, these values 
were used to estimate the parameters of a beta distribution 
[Ref. 17, p. 196] . 

Observations of the daily sewage generation rate per 
man were then generated from this beta distribution for each 
day of the simulated week. Once a daily average was gen- 
erated it was converted to an hourly average for a ship 
with 240 men. (The number of men was arbitrary and could 
be read in by the programmer.) Finally, it was translated 
to the empirical distribution by conversion ratios of the 
hourly relative value to the relative mean. For example 
at 1200, if the beta realization was 15 gal/man/day then: 
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Relative Value 

Hourly Realization = — 

Relative Mean 

Daily Beta Variate 
X — X No. of Men 
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11.3 15 

Realization = x — x 240 = 232 gal 

7.3 24 



( 12 ) 



Such values were the hourly realizations of X(t) used 
as the input to the basic model described in subsection 1. 
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IV. RESULTS 



The results of simulations one and two are presented on 
pages 32 through 35. Figure 5a and 6a are graphs of the 
average maximum level of sewage held as a function of pro- 
cessing rate for the five policies tested by simulations one 
and two respectively. Figures 5b and 6b represent the 
average total time the processor must r\an in a week's period 
(of a possible 168 hour maximum) , again as a function of 
processing rate for each of the five policies. The results 
agree well with intuition. 

Looking first at the results common to both simulations, 
one observes that the maximiam average sewage level has an 
inverse relationship to the processing rate. The percent 
savings in tank size decreases with higher processing rates, 
and, for rates greater than 210 gal/hr, the savings appear 
to be negligible. For policies in which processing time 
can vary, the total weekly processing time is also inversely 
proportional to the processing rate, and, over the range of 
rates considered, the function is nearly linear. 

The different policies considered have a significant 
effect on sewage level and processing time. Policy two 
describes a limit on minimixm sewage level obtainable for 
the rates considered. The curve of maximum sewage level 
increases with each higher policy number up to policy five. 
Conversely, for average processing time the trend is just 
the opposite as was expected. Policy two requires the 
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maximum processing time and each succeeding policy number 
requires less. For the policy of strictly holding sewage 
(policy 1) , the average sewage accumulated in a week's 
period (asstiming the ship to be a destroyer with a comple- 
ment of 240 men) was 21,100 gallons, with a variance of 
8.2 X lO'' gallons^ for simulation one, and 21,080 gallons 
with a variance of .23.4 x 10® gallons^ for simulation two. 

The average values for a 12 hour holding period were 3,057 
gallons for simulation one and 2,700 gallons for simulation 
two. 

Certain additional observations could be made that were 
peculiar only to simulation one. First, a definite change 
in the slope of the curve for maximum tank level vs pro- 
cessing rate occurs at 150 gal/hr, with little additional 
savings in tank size being realized beyond that. For pro- 
cessing rates above 150 gal/hr, the average maximum tank 
level was never above 2,500 gallons. The variance/mean 
ratio for maximtom tank level is on the order of from 2 to 10. 
It is usually highest for a processing rate of 140 gal/hr, 
and drops off immediately thereafter. For policy four the 
ratio is much lower than for other policies. For processing 
time vs processing rate the variance/mean ratio is on the 
order of 10“^ , and policies 3 and 4 show little appreciable 
difference in time savings. 

For simulation two, it can be seen that, with policy 
two as a reference, the curve of maximum level vs processing 
rate levels off more rapidly for policies 3 and 4 beginning 
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at approximately 160 gallon/hr. The curves for policies 3 
and 4 follow the same contour with policy 3 having a maxi- 
mxim sewage level on the order of 1000 gallons less than 
policy four for all rates. The variance/mean ratio for 
maximum sewage level is on the order of 10^ to 10^ at 140 
gal/hr and decreases as processing rate increases. At 
higher .processing .rates ..the ratio is significantly less for 
policy four (approx. 20) than for other policies. For pro- 
cessing time vs processing rate, policies three and four 
have the same slope, with policy four showing only a slight 
savings in time of about 5 hrs for the week. The variance/ 
mean ratio for either policy is again very small , on the 
order of 0.5. 

A comparison of the two simulations reveals that, al- 
though the average results of simulation one are on the 

» 

same order as those of simulation two, the variance of the 
results is much smaller for the former. This is not supris- 
ing, since simulation one is based on the Poisson process 
in which the variance of the asymptoptic distribution of 
N(t) will equal the mean. In the case of a non-homogeneous 
process in which the parameters do not differ greatly, it 
can also be expected that as the number of replicates in- 
creases the sample variance will approach some value on the 
same order of magnitude as the mean. For simulation two, 
the realizations of the daily average sewage rate per man 
can vary from 11.3 gal/man/day to 20.0 gal/man/day which is 
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then imposed on an empirical distribution with even greater 
variance. The result is a process with reasonably large 
variance; certainly much greater than that resulting from 
simulation one. 
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Figure 5a. Maximum Tank Level vs Processing 
Rate for Simulation One. 
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Processing Rate (Gal/Hr) 
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Processing Rate (Gal/Hr) 



Figure 6a. Maximum Tank Level vs 
Processing Rate for 
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Processing Rate (Gal/Hr) 



V. CONCLUSIONS AND RECOMMENDATIONS 

A. INTERPRETATION OF RESULTS 

In order to prescribe the most feasible combination of 
processing rates, policies and holding tank size, more 
information regarding the cost fxinction of each particular 
parameter would have to be developed. If costs were as- 
sumed' to follow some 'increasing linear function, of pro- 
cessing rate, holding tank size and processing time, then 
the policy and rate which would minimize a measure of ef- 
fectiveness taking all such costs into account would be 
expected to be most feasible. 

Based on the results of Chapter IV and with the under- 
standing that cost functions must be applied, it appears 
that a system consisting of a holding tank of 2,700 gallon 
capacity for sewage combined with a processing unit capable 
of handling 160 to 170 gallons of sewage per hour would be 
most beneficial. Further such a system should be operated 
under a processing policy similar to that of policy 3. 

It was contended in Chapter I that the Navy's design 
parameter of 30 gallons/man/day was over estimated. The 
simulations, which were developed using the best data that 
could be obtained, bear out this conclusion. More specifi- 
cally, if this current value were used to construct a 12 
hour holding tank for a destroyer of 240 men, the minimum 
tank size would be 3600 gallons. The values obtained 
through simulation, which were measured during the period 
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when traffic flow was the greatest (0500 to 1700) , were, 
again, 3,050 gallons for simulation one, and 2,700 gallons 
for simulation two. Even with the addition of a 200 gallon 
safety factor, it appears evident that a savings of between 
400 to 700 gallons could be realized for each 12 hour 
holding tank installed. Even greater savings would be 
expected..for larger tanks. 

The capacity of a holding tank based on the quantity of 
diluent (flushing liquid) used can also be computed using 
an equation presented in a report by the Maritime Research 
Information Service [Ref. 13]. Namely: 

C = 1.25 X M X (B+D) (13) 

where 

C = tank capacity 

M = man days of storage 

B = daily per capita body waste = 0.437 gallons 

D = diluent in gallons/man/day . 

If the value of diluent is taken to be 15 gal/man/day, 
the unmodified mean of the data obtained by MEL Report 
346/64, then the recommended holding tank size would be 
2,320 gallons. The values obtained by simulation would 
thus appear conservative; however, it is important to 
realize that they also account for peak arrival periods. 

B. RECOMMENDATIONS 

The simulations described in this thesis offer sufficient 
merit for use in the design of sewage treatment systems. It 
is recommended that a study be made to obtain further 
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statistics on the arrival pattern of head facilities aboard 
a Naval ship. Such new data could then replace the empirical 
distribution developed for simulation two, or be used to 
revise the intensity function of simulation one. This 
should lead to even more meaningful results than those 
presented . 

It is further recommended that design parameter for 
daily per capita sewage generation be reviewed. If it is 
decided that it is in fact overestimated, then it is recom- 
mended that 12 hour holding tanks of the size presented in 
this report be tested on board destroyer-sized ships, and 
that the simulations be employed to develop holding tanks 
for larger ships as well. 
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APPENDIX A 



CALCULATION OF ARRIVAL PARAMETERS 

As pointed out in the body of the thesis , the mean value 
function was used to estimate the best values of Xj^ (i= 
1,2,3). This was done by computing average daily and weekly 
sewage generation first by using the average values pre- 
sented in MEL Report 346/64 [Ref. 19] , and second by use of 
policy one and simulation two. The values obtained from 
weekly observations of sewage build-up as estimated by 
simulation two were used since 1 ) they were considered more 
representive and 2 ) they would also serve as a means of 
normalizing both simulations so that over long periods of 
time the two simulations would give similar results for 
strict holding. Average daily and weekly sewage generation 
was then divided by 4.9 gal/arrival to determine the number 
of arrivals/period. 

The results were : 

Weekly average = 21,000 gal/week/4.9 gal/arrival 

= 4260 arrivals/week 

similarly : 

daily average = 730 arrivals/day . 

These were then equated to the mean value function which 
was 



M(t) = 5(m(t)) + 2 (m(t')) = 5 [ X j ( t j ) +X 2 ( 1 2 -t i ) +X 3 { t-tj ) ] 
+ 2[X2(t/)+X3(t-t/)] (A-1) 
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for the weekly value , and 



m(t) = (tj)+X 2 (t 2 ’-ti)+X 3 (t-t^) (A-2) 

for the daily value. 

Substituting in the values for tj^ where t = 24 hours = 
288 five minute increments and setting them equal to the 
weekly and daily sewage build-ups, one gets: 



weekly : 


66 OX 1 + 768 X 2 + 588 X 3 


= 4260 


(A-3) 


daily : 


132Xi + 72 X 2 + 84 X 3 = 


730 


(A-4) 


where : 


tj = 132, t 2 = 204, 


t = 288, t^' = 132. 





In order to solve the above simultaneous equations a 
third equation was needed. This was obtained from weekends 
and computing a third average build-up for that period to 
obtain : 

2 m(t') = 2,100/4.9 = 420 arrivals/weekend 
or 

264X2 312 X 3 = 420. (A-5) 

Solving (A-3) through (A-5) simultaneously one obtains 
the final estimates of X^ (i = 1,2,3) namely: 

A 

Xj = 4.5 arrivals/5 min 
X 2 = 1.0 arrivals/5 min 
X 3 = .5 arrivals/5 min 

which agree well with intuition. 
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APPENDIX B 



METHOD OF ESTIMATION FOR AN EMPIRICAL DISTRIBUTION 

In order to obtain the hydraulic loading pattern of 
water closets and urinals aboard Naval ships, it was nec- 
essary to resort to a subjective estimation technique. 

The technique used is the same as that for obtaining job 
time estimates for PERT scheduling computations [Ref. 14] . 
Such a technique is employed when one is unable to do any 
statistical sampling. 

It was, therefore, decided to survey a cross-section 
of experienced Naval personnel, both officers and enlis^ted, 
for their subjective estimates of the relative usage of 
head facilities for each hour of a 24 hour period. Their 
estimates were made on a relative scale of 0 to 20, (0 in- 

dicating the facility was not used at all, 20 indicating 
it was used to capacity) . Each person was to give a high, 
most likely, and low estimate, independent of estimates 
made by others . 

It has been "historically accepted" that job duration 
times described by such estimates are beta distributed over 
a finite interval given by the high and low estimate. The 
same distribtuion can be assumed to apply to estimates of 
the arrival pattern. The expected value of such a distri- 
bution can then be approximated by: 

Ug = (a+4m+b) /6 (B-1) 
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and the variance by; 



= I(b-a)/6j2 (B-2) 

where ; 

a = low (optimistic) estimate 
m = most likely (mode) estimate 
b = high (pessimistic) estimate. 

Using the expected value, , as the best estimate of each 
individual. All such values were then averaged to give the 
overall best estimate for each hour of the 24 hour period. 
The results are the empirical distribution plotted in the 
body of the thesis (Figure 4a) . 
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APPENDIX C 



ESTIMATION OF BETA PARAMETERS 



A beta distribution that is transformed to some interval 
other than (0,1) is given by: 



r (A+B) 



(b-a)l-A-B (,p_^)A-l (b-T)B-l^ b>T>a 



f(T) = 



r(A) r(B) 

0 elsewhere 

where T is the transformed beta variable given by: 



(C-1) 



T = a+(b-a)X (C-2) 

X being a (0,1) beta variable and a,b being the lower and 
upper bound of the transformed variable. If the high (b) , 
low (a) and most likely (m) values of the distribution can 
be determined, then the transformed beta parameters can be 
related to (0,1) beta parameters according to: 

A = 46E(X) IE(X)-E(X)^ - (C-3) 

A-AE(X) 

B = (C-4) 

E(X) 

where 



4m+b-5a 

E(X) = (C-5) 

6 (b-a) 

Using the values obtained from MEL Report 346/64 i.e.: 



a = 11.3 gal/man/day 
m = 14.8 gal/man/day 
b = 20.0 gal/man day. 
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Then the values for the parameters can be solved for, using 
the above equations to obtain: 

A = 4.3 B = 5.7 

These values are then used to generate a beta variable, X, 
by first generating two independent gamma variates yj~G(l,A) 
and ¥ 2 ( 1 , 6 ) and computing the ratio 



X = 


(C- 6 ) 


It then remains 


only to transform X to a realization of T 



using equation (C-2) . 
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APPENDIX D 



USE OF SIMULATION ONE AND TWO 



The two simulations were programmed in Fortran IV for 
use on the IBM 360/67 computer. Both simulations allow the 
option of plotting the graph of the sewage build-up (in- 
ventory) for each hour of the week. (Examples of such 
plots are given under computer output.) If one uses the 
option, however, he must insure that the appropriate JCL 
cards are used for his particular hardware. 

Only one data card is needed for either simulation. 

All that must be provided are the parameters to be tested, 
an indication of print routine to be used and a flag for 
calling the plot routine. For each simulation the format 
of the data card is as follows : 

a) Simulation One. The following variables must 



be supplied: 



Variable 



Cols on Card 



Description 



Plan 1 



Plan 2 



6-10 



0-5 



As described in the body 
of the thesis, these are 
the arrival rates/5 min 
that are used for the 



Plan 3 



11-15 



Poisson generator. 



Alpha 1 



16-18 



The initial processing 
rate seed, in gal/5 min 
processed. It is in- 
creased by 10 gal/hr as 
soon as it is read in . 



NPLAN 



19-20 



A two digit integer 
specifying the processing 
plan to be tested. 
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NOPRIN 



21-22 



NO DRAW 23-24 



ZSTOP,ZSTART 25-20 

30-34 



A two digit integer to 
indicate which of four 
print routines is to be 
used. (See below for 
discription of these 
routines . ) 

A flag to indicate if 
plot routine is desired 
(any 2 digit number) or 
not (use zero) . 

These are the level 
parameters for policy 4 
in F5 . 1 format . 



TSTOP 



34-36 A 3 digit integer repre- 

senting the stopping 
time for policy 5. 



For simulation two the same values are read in with the 



exception of the Poisson parameters, instead the number of 
men representing the size of ship is read in as a 3 digit 
integer . 

For the print routines it should be noted that if full 



printout is used that this can add considerably to the run- 



ning time of the program. 



Depending on the read in for 



NOPRIN , the routines are : 
NOPRIN 
1 



2 



3 



ROUTINE 



Print only the final statis- 
tics for ZMAX matrix and the 
processing time matrix. 

Print above plus the times 
that the processor stop and 
restarted during the week. 
Print also the processing 
rate , run nvimber and maximum 
tank level on each run. 

Print the above plus the 
entire Z(t) matrix (also the 
X(t) matrix for simulation 2). 
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4 Print the above plus the 

X(t) matrix and arrivals 
matrix. (Simulation One) 

In conclusion, it should be pointed out that simulation 
one requires approximately 2000 random numbers to be gen- 
erated per replication; it requires considerably more com- 
puter time than simulation two. This is inpart compensated 
for, however, by its smaller variance, so that fewer runs 
are needed. 
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K-SCRLE-2.D0E+Di UNITS INCH. 
y-SCnLE---2. 00E+D2 UNITS INCH. 

TRNKSIHi 



COMPUTER PROGRAM 



* TANK SIMULATION ONE * 

OU ^iU ^1# «JU WU «JU «JU» OU WU OU «JU %V «4^ OU «4U «JU ^4^ «JU «4U ^U OU ^L# «1# «4U ^U «JU «JU %A^ ^L# «JU «JLr 

THIS IS THE FIRST OF TWO SIMULATIONS DEVELOPED 3Y J.O.MINER» 
JR. TO DESCRIBE THE PROCESS OF GENERATING AND TREATING 
SEWAGE ON BOARD NAVAL SHIPS. IT EMPLOYS A NON-HOMOGENEOUS 
POISSON PROCESS TO DESRIBE THE ARRIVAL OF SEWAGE TO A HOLD- 
ING TANK OVER THE COURSE OF THE DAY. IT THEN TESTS FIVE PRO- 
CESSING POLICIES FOR ANY PROCESSING RATE DESIRED. THE SIZE 
OF THE SHIP IS ARBITRARY AND IS READ IN BY THE PROGRAMMER. 

» >1^ 4#^ ^1^ 4i^ >1^ #1^ >1^ 4^lf^ >1^ ^Y^ 4i^ «Y^ ^1^ 

jjt * 

-JO -O 

THE FOLLOWING STATEMENTS INITIALIZE THE PROGRAM AND READ In" 
THE VARIOUS PARAMETERS TO BE TESTED, INCLUDING THE PARAMETERS 
FOR THE POISSON DISTRIBUTION 

COMMON I,J,ITIME,NOPRIN 

DIMENSION Z(7,288) ,X(7,288 ) ,NUM(7,288) 

DIMENSION XORD( 168) ,YORD( 168), ZMAXd 0,20) , TPROSS ( 10 , 20 
REAL=!'8 TITLE(12)/'TANKSIM1 • , 11=!'' •/ 

REAL LABEL/* •/ 

INTEGER TSTOP 

READI5 ,100) PLAMl, PLAM2,PLAM3, ALPHA1,NPLAN,NQPRIN, 
CNODRAW,ZSTOP,Z START, TSTOP 
100 F0RMAT(3F5.3,F3.1,3I2,2F5.1,I3) 

IX=4321571 

RUNS ARE MADE FOR NINE PROCESSING RATE IN INCREMENTS OF 10 
GALLCN/HOUR 

DO 6000 M=l,9 
ALPHAl=ALPHAl+.833 
DO 5000 L=l,10 
IFINOPRIN.lt. 2) GO TO 155 
WRITE(6,150)L 

150 FORMAT {• 1* ,20X, • RUN NO. *,I2,//) 

VARIOUS ARRAYS ARE ZEROED FOR EACH RUN 

^ 4= 

155 ZMAXI M,L)=0.0 
ALPHA=0.0 
ITIME=0 
DO 30 1=1,7 
DO 20 J=l,288 
Z(I,J)=0.0 
X( I , J) = 0.0 
NUM{ I , J) = 0 
20 CONTINUE 
30 CONTINUE 

IX = IX + 342 

3{c ^ 

A POISSON ARRIVAL IS THEN GENERATED FOR EACH FIVE MINUTE 
INCREMENT OF THE SIMULATED WEEK. DIFFERENT PARAMETERS ARE 
USED ACCORDING TO THE TIME OF DAY OR WEEK. 

DO 600 1=1,5 

PLAM = PLAMl 

DO 500 K=l,132 

CALL P01SSN(PLAM,IX,N) 

500 NUM(I,K)=N 
PLAM =PLAM2 
DO 510 K=133,204 
CALL POISSNIPLAM, IX,N) 

510 NUM{ I ,K)=N 
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PLAM = PLAM3 

DO 520 K=205,288 

CALL P01SSN(PLAM, IX,N) 

520 NUM(I,K)=N 
600 CONTINUE 

DO 750 1=6,7 

PLAM= PLAM2 

DO 700 K=l,204 

CALL P01SSN(PLAM,IX,N) 

700 NUM( I ,K)=N 
PLAM=PLAM3 
DO 720 K=205,288 
CALL P01SSN(PLAM, IX,N) 

720 NUM(I,K)=N 
750 CONTINUE 

DO 1500 1=1,7 

THE ARRIVALS ARE THEN CONVERTED TO A QUANTITY OF SEWAGE ANd" 

The sewage processed in accordance with one 0F-‘THE five pol- 
icies. 

if. 

IF( (NPLAN.EQ.l ) .OR. (NPLAN.EQ.4) ) GO TO 90 
ALPHA=ALPHA1 
90 DO 1000 J=l,288 

X(I ,J)=NUM(I,J)=!^4.94 

IF( ( J. EQ. 1 I.AND. (I .EQ.l ) )G0 TO 855 

IF( ( J . EQ. 1) .AND. ( I .GE.2) )G0 TO 860 

Z(I,J)= Z( I,J-1)+X( I, J )-ALPHA 

GO TO 880 

855 IF (NPLAN.EQ.4) ALPHA=ALPHA1 
Z( I, J)=X( I,J)-ALPHA 
GO TO 880 

860 Z(I,J)= Z(I-1 ,288)+X( I , J)-ALPHA 
830 IF(Z(I,J).LT.0.0)Z(I,J)=0.0 

IF(ZMAX(M,L) .LT.Z( I,J) ) ZM AX ( M , L ) = Z ( I , J ) 

TANK=Z(I,J) 

IF(NPLAN.EQ.l) GO TO 1000 
IF(NPLAN.EQ.2) GO TO 1000 

IF (NPLAN.EQ.3)CALL PL AN3 ( T ANK , AL PH A 1 , AL PH A , 8900 ) 

IF (NPLAN.EQ.4) CALL PLAN4(TANK, AL PH A 1 , ALP HA , ZSTOP , 
CZSTART , 8900) 

IF(NPLAN.EQ.5) CALL PLAN 5( ALPHA 1 , ALPHA , TSTOP , 8900 ) 

900 CONTINUE 
1000 CONTINUE 
1500 CONTINUE 

THE WEEKLY PROCESSING TIME IS THEN COMPUTED 

IF( NPLAN.EQ.l) ITIME=2016 
IF(NPLAN.EQ.2) ITIME=0 
IF(NPLAN.EQ.5) I T I ME= ( 288-TSTOP ) «7 
TPROSS(M,L) =168. -I TIME/12. 

IF DESIRED, THE WEEKLY SEWAGE INVENTORY MAY BE PLOTTED. 

IF(NODRAW.EQ.O) GO TO 2900 

INDEX=0 

DO 2000 1=1,7 

DO 1900 J=12,288,12 

INDEX=INDEX+1 

YORD( INDEX)=Z( I, J) 

1900 CONTINUE 
2000 CONTINUE 

DO 2100 J=l,168 
2100 XORD( J )=J 

CALL DRAW(168,X0RD,Y0RD,0, 0, LABEL, TITLE, 20.0,300.0,0, 0 
ILAST) 

REALIZATION OF EITHER SEWAGE OR ARRIVALS FOR EACH FIVE MIN- 
UTE INTERVAL MAY THEN BE PRINTED AS DESIRED. Z VALUES ARE 
CUMMULATIVE HOURLY 
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2900 IF(N0PRIN.LT.2) GO' TO 5000 

WRITE(6,3500)2MAX(M,L) , ALP HA 1 , TP RO SS ( M , L ) 

3500 FORMATCO' ,20X,'ZMAX=' ,F8.1,10Xt 'PROCESSING RATE = 
CF6.1 ,10X,' PROCESSING TIME = '»F7.2,//) 

2950 IF (NOPRIN.lt. 3) GO TO 5000 

WRITE (6, 3000) ((Z(I»J),I=l,7)tJ=12,288,12) 

3000 F0RMAT('0 Z MATRI X ' // ( 1 X, 7 ( F 8 . 1 ,4X) ) ) 

4000 IF(N0PRIN.LT.4) GO TO 5000 

WRITE (6, 3100) ( (X( I ,J) , 1 = 1,7 ) ,J = 1,288 ) 

3100 FORMATCl X MATRI X • // ( IX , 7 ( F 5. 2 ,3X) ) ) 

WRITE(6,3200) ( (NUM( I, J ) , 1=1,7) , J=1 ,288) 

3200 FORMATCl ARRIVALS MATR IX • // ( IX ,7 ( 12, 3X ) ) ) 

5000 CONTINUE 
6000 CONTINUE 

STATISTICS ARE THEN COMPUTED FOR MAXIMUM SEWAGE LEVEL AND 
FOR PROCESSING TIME. THE FINAL RESULTS ARE THEN PRINTED 

DIMENSION ‘ZSUM(20) ,ZAVE(20) ,ZVSUM(20) ,TSUM(20) ,TAVE(20 
C) ,TVSUM(20) ,ZVAR(20) ,TVAR( 20) 

DO 7200 J=1,M 
ZSUM( J )=0. 

ZVSUM( J)=0. 

TSUM( J)=0. 

TVSUM( J )=0. 

DO 7000 1=1, L 

ZSUM( J)=ZSUM( J)+ZMAX( J ,I ) 

TSUM(J)=TSUM(J)+TPROSS(J,I ) 

7000 CONTINUE 

ZAVE(J)=ZSUM( J) /L 
TAVE( J )=TSUM(J )/L 
DO 7100 1=1, L 

ZVSUM( J)=ZVSUM( J)+( (ZMAX( J,I)-ZAVE(J) )**2) 

TVSUM( J)=TVSUM(J)+((TPROSS(J,I )-TAVE(J) )^^Z) 

ZVAR( J) =ZVSUM( J)/L 
TVAR( J )=TVSUM( J) /L 
7100 CONTINUE 
7200 CONTINUE 

WRITE(6,8000) NPLAN 

8000 FORMATCl ZMAX MATRIX FOR POLICY NO. ',12//) 

DO 8200 J=1,L 

WRITE(6, 8100) ( ZMAX( I ,J ) , 1 = 1 ,M) 

8100 FORMAT (IX, 10( F8 .2,3X) ) 

8200 CONTINUE 

WRITE (6, 8 500) (ZAVE( J) , J=1,M) 

8500 FORMATCO AVERAGE ZMAX ' / / ( 3X , 10 ( F8 . 2 , 2X ) ) ) 

WRITE! 6,8600) ( ZVAR( J) , J = 1,M) 

8600 FORMATCO ZMAX VARIANCE ' // ( 3X, 10 ( F 10. 2 , 1 X ) ) ) 

WRITE! 6 ,9000) NPLAN 

9000 FORMATCl PROCESSING TIME MATRIX FOR POLICY NO. ', 
CI2// ) 

DO 9200 J=1,L 

WRITE! 6,8150 ) (TPROSS! I ,J) ,I=1,M) 

8150 FORMAT!1X,10!F7.2,3X)) 

9200 CONTINUE 

WRITE (6, 9300) (TAVE! J) , J=1,M) 

9300 FORMATCO AVERAGE T IME ' / / ! 3X , 10 ! F 6. 2 , 2X ) ) ) 

WRITE! 6,9400) (TVAR(J) , J=1,M) 

9400 FORMATCO TIME VARIANCE ' // ! 3X, 1 0 ! F 1 0. 2 , 1 X) ) ) 

STOP 

END 

•JU -JU 

THE FOLLOWING SUBROUTINE GENERATES POISSON VARIATES FDR EACH 
FIVE MINUTE INCREMENT. A TABLE LOCK-UP METHOD WAS USED 
* ❖ 
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SUBROUTINE POISSNI PLAN , I X,N ) 
COMMON I ,J , ITIMEtNOPRIN 
INTEGER N 
A=EXP(-PLAM) 

B=A*PLAM+A 
C=(At>PLAM*=«’2)/2+B 
D=( A*PLAM**3) /6+C 
E=( A=!=PLAM=l'=!=4)/24+D 
F=( A*PLAM**5) /120+E 
G=( A’!'PLAM*=!=6)/720+F 
H=( A*PLAM**7)/5040.+G 
P=( A*PLAM**8) /40320.+H 

5 CALL RANDU( IX, IY,YFL ) 

IX=IY 

R= YFL 

IF(R.LE.A) GO TO 6 
IF (R.LE.B) GO TO 7 

IF (R.LE.C) GO TO 8 

IF (R.LE.D) GO TO 9 

IF (R.LE.E) GO TO 10 

IF (R.LE.F) GO TO 11 

IF (R.LE.G) GO TO 12 

IF (R.LE.H) GO TO 13 

IF (R.LE.P) GO TO 15 

N=9 

GO TO 14 

6 N=0 

GO TO 14 

7 N=1 

GO TO 14 

8 N=2 

GO TO 14 

9 N=3 

GO TO 14 

10 N=4 

GO TO 14 

11 N=5 

GO TO 14 

12 N=6 

GO TO 14 

13 N=7 

GO TO 14 
15 N=8 

14 RETURN 
END 



SUBROUTINE PLAN3 (T ANK , ALFA 1 , RATE ,* ) 

THIS POLICY PROCESSES SEWAGE UNTIL 1600 AFTER WHICH TIME IT 
THEN EXAMINES THE LEVEL OF THE TANK AND IF BELOW 50 GAL. IT 
STOPS PROCESSING UNTIL 0500 THE FOLLOWING DAY. ON WEEKENDS 
PROCESSING STOPS ANYTIME THE LEVEL IN THE TANK REACHES 50 GA 

COMMON I ,J ,ITI MEtNOPRIN , 

IF(( I.EQ.lJ.AND.(J.EQ.l) ) IFLAG=0 ' 

IF(N, OPRIN.lt. 2) GO TO 6 

IF({ I.EQ.D.AND. (J.EQ.l) ) WRITE (6, 5) 

5 FORMATCO PROCESSING POLICY NO. 3 WAS USED ON THIS 
CRUN' ,//) 

6 IF((J.EQ.1).AND.(IFLAG.EQ.D) I TIME=I TI ME+288-JST0P 
IF(J.EQ.l) IFLAG=0 

IF( I .LE.5) GO TO 10 
IF(TANK.LE.50. ) GO TO 30 

IF (( RATE. LT. 0.01 ) .AND. (TANK. GE .500. ) ) GO TO 20 
GO TO 40 

10 IF(J.LE.132) GO TO 40 
IF(TANK.LE.50. ) GO TO 30 
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GO TO 40 
20 RATE=ALFA1 

IF( IFLAG.EQ.O) GO TO 40 
ITIME=ITIME+J-JSTOP 
IF(NCPRIN.LT.2) GO TO 26 
WRITE(6,25)I ,J 

25 FORMATC CONTINUED PROCESSING AT T=' , II , • , 13 , // ) 

26 IFLAG=0 
GO TO 40 

30 RATE=0.0 

IFdFLAG.EQ.DGO TO 40 
JSTOP=J 

IF(N0PRIN.LT.2) GO TO 36 
WRITE(6,35)I ,J 

35 FORMATI'O STOPPED PROCESSING AT T= S 1 1 , ' - • t I 3 , / / ) 

36 IFLAG=1 

40 IF( n.EQ.7).AND.(J.EQ.288) ) GO TO 45 
GO TO ,50 

45 TF{ IFLAG.EQ. 1) IT IME = ITIMB+268-' JSTOP 
50 RETURNl 
END 



SUBROUTINE PLAN4 ( TANK , ALFAl , RATE , STOP , START,*) 

THIS IS A LEVEL CONTROL POLICY WHICH CONTINUOUSLY EXAMINES 
T&& LSOEL OF SEWAGE IN THE TANK. IF BELOW SOME SET LEVEL 
(ZSTOP) THEN PROCESSING STOPS UNTIL THE SEWAGE LEVEL REACHES 
SOME UPPER VALUE (ZSTART) AT WHICH TIME PROCESSING CONTINUE 
* 

COMMON I ,J,ITIME,NOPRIN 

IF(( I.EQ.D.AND. (J.EQ.l) ) IFLAG = 0 

IF(N0PRIN.LT.2) GO TO 6 

IF(( I.EQ.D.AND. (J.EQ.l) ) WRITE (6, 5) 

5 FORMATCO PROCESSING POLICY NO. 4 WAS USED ON THIS 
CRUN' ,//) 

6 IF( ( J. EQ. 1). AND. (IFLAG.EQ. 1) ) I T IME=I TI ME+288-JSTOP 
IF(J.EQ.l) JSTOP=0 

IF(TANK.GT.STOP) GO TO 20 
RATE=0.0 

IFdFLAG.EQ.DGO. TO 40 
JSTOP=J 

IF(N0PRIN.LT.2) GO TO 11 
WRITE(6,10) I ,J 

10 FORMATI'O STOPPED PROCESSING AT T= • , II , , 13, // ) 

11 IFLAG=1 
GO TO 40 

20 IFITANK.LE. START) GO TO 40 
RATE=ALFA1 

IF ( I FLAG.EQ.O) GO TO 40 
ITIME=ITIME+J-JSTOP 
IF(N0PRIN.LT.2) GO TO 31 
WRITE! 6,30) I ,J 

30 FORMATC CONTINUED PROCESSING AT T=* , 1 1 , , I 3 , // ) 

31 IFLAG=0 

40 IF( (I . EQ.7) .AND. (J.EQ.288) ) GO TO 45 
GO TO 50 

45 IFdFLAG.EQ.l ) I T I ME= IT I ME + 288- J ST OP 
50 RETURNl 
END 



SUBROUTINE PL AN5 ( ALFAl , RAT E ,T I ME , * ) 

THIS IS A TIME CONTROL POLICY. IT ALLOWS PROCESSING FOR A 
FIXED PORTION OF THE DAY ONLY 

COMMON I , J , ITIME,NOPRIN 

INTEGER TIME 

IF(N0PRIN.LT.2) GO TO 6 

IFd I .EQ.l ) .AND. ( J .EQ. 1) ) WRITE! 6, 5) 

5 FORMATI'O PROCESSING POLICY NO. 5 WAS JSEO ON THIS 
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CRUN' ,//) 

6 IF(J .LE.TIME) GO TO 10 
RATE=0. 0 
10 RETURNl 
END 



SUBROUTINE RANDU ( I X, lY , YFL) 

THIS ROUTINE GENERATES RANDOM NUMBERS FOR USE IN THE MAIN 
PROGRAM 

COMMON I ,J, ITIMEtNOPRIN 

IY=.IX^65539 

IF(IY)5,6,6 

5 IY=I Y+2147483647+1 

6 YFL=IY 

YFL=YFL* .46566 13 E-9 

RETURN 

END 



Note: Comment cards denoted by and including * should be 

removed prior to running programs. 
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* TANK SIMULATION TWO * 

THIS SIMULATION WAS DEVELOPED BY JOHN 0 MINER»JR TO DESCRIBE 
THE PROCESS OF GENERATING AND PROCESSING SEWAGE ABOARD A 
NAVAL SHIP. IT ASSUMES THAT THE DAILY AVERAGE SEWAGE GENERA- 
TION RATE PER CAPITA IS BETA DISTRIBUTED OVER THE INTERVAL 
11.3 GAL/MAN/DAY TO 20.0 GAL/MAN/DAY. IT GENERATES DAILY 
OBSERVATIONS OF THE DAILY RATE AND CONVERTS THEM TO HOURLY 
VARIATES .BASED ON AM EMPIRICAL ARRIVAL DISTRIBUTION. THE 
HOURLY OBS'ERVATI'ONS THEN SERVE AS 'THE 'INPUT TO A SEWAGE 
HOLDING TANK FROM WHICH THE SEWAGE IS PROCESSED ACCORDING TO 
ONE OF FIVE SAMPLE PROCESSING POLICIES WHICH ARE ARE REPRE- 
SENTED BY INDIVIDUAL SUBROUTINES. OVER A WEEK'S PERIOD OF 
SIMULATED TIME, THE MAXIMUM SEWAGE LEVEL OBSERVED FOR UP TO 
100 REPLICATES IS RECORDED ALONG WITH THE TOTAL PROCESSING 
TIME. STATISTICS ARE THEN COMPUTED AND THE THE RESULTS PRINT 
ED ACCORDING TO ONE OF FOUR POSSIBLE PRINT ROUTINES. 

^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ■JU ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 

^ ^ ^ ^ ^ ^ ^ ^ ^ 41^ ^ ^ Sf. ^ Ac- ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 'C’ ^ "T‘ ^ ^ 'f‘ 'C’ 

^ 2*2 

THE FOLLOWING STATEMENTS INITIALZE THE PROGRAM. THE DATA 
STATEMENT READS IN THE CONVERSION FACTORS FOR THE EMPIRICAL 
DISTRIBUTION 

^ 2{C 

COMMON I X, I , J , ITIME,NOPRIN 
DIMENSION C0NVER(24) ,X(7,24),Z(7.24) 

DIMENSION DAILY (7) , HOURLY! 7) ,Y(7),Y1(7) 

DIMENSION X0RD(168) ,Y0RD(168) 

DIMENSION ZMAX( 10, 100) , IPROSS! 10,100) 

REAL*8 TITLE(12)/'TANKSIM2' , 11*' ’/ 

REAL LABEL/' •/ 

DATA CONVER/. 756, 1.9 7, 1.82, .963, .812, 1.09, 1.80, 1.61, .9 
D63,. 82 5, 1.47, 1.69, 1.31, 1.04, .839,. 77 0, .880, .770,. 687, 
E.495, .344, . 234, . 440, . 413/ 

THE PARAMETERS TO BE TESTED ARE THEN READ IN 

READ(5,20)MEN, NP LAN, ALPHA2,ZST0P ,Z START, NOPR IN, NODRAW 
C,KSTOP 

20 F0RMAT(I3,I2,F4. 1,F5.2,F5. 1,312) 

35c ;{c 

THE FIVE POLICIES ARE THEN TESTED SEQUENTIALLY FOR 100 REP- 
LICATIONS EACH. THE SAME RANDOM NUMBER SEED IS USED FOR EACH 
POLICY 

DO 9999 NRUN=1,4 
IX=21111111 
NPLAN=NPLAN+1 
ALPHA1=ALPHA2 

3{C SJC 

EACH POLICY IS TESTED FOR NINE RATES IN INCREMENTS OF 10 GAL 
245 ^ 

DO 6000 M=l,9 
ALPHA1=ALPHA1+10.0 
DO 5000 L=l,100 
IX=IX+112 

IF(NOPRIN. LT.2) GO TO 155 
WRIFE(6, 150)L 

150 FORMAT (' 1' ,20X, 'RUN NO. ',12,//) 

155 ZMAX(M,L)=0.0 
ALPHA=0.0 
ITIME=0 
DO 30 1=1,7 
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DO 25 J=l,24 
Z(I,J)=0.0 
X( I ,J) = 0.0 
25 CONTINUE 
30 CONTINUE 

SEVEN BETA VARIATES ARE GENERATED FOR EACH DAY OF THE WEEK 

DO 50 1=1,7 
Y( I)=GAM(4.3,1.,1. ) 

YKI )=GAM(5.7,1.,1.) 

BETA2=Y( I ) /(Y{ I )+Yl( I ) ) 

DAILY! I )=BETA2*8.7+11.3 
50 CONTINUE 

IF(N0PRIN.LT.2) GO TO 56 
WRIT E( 6, 55) (DAILY! I ) ,1=1,7) 

55 FORMATCO DAILY BETA VARIATES: • // ( IX , 7 ! F5 .2 , IX ) ) ) 

THE DAILY VARIATES ARE CONVERTED TO HOURLY VALUES AND 
THEN TRANSFORMED TO THE EMPIRICAL DISTRIBUTION 

' 56 DO 60 1=1,5 

60 H0URLY!I)=!DAILY!I)*MEN)/24 
DO 65 1=6,7 

65 HOURLY ! I ) =!DAILY! I )*MEN/4) /24 
DO 80 1=1,7 
DO 70 J=l,24 

X! I , J) =CONVER! J)*HOURLY! I) 

70 CONTINUE 
80 CONTINUE 

DO 1500 1=1,7 

IF! (NPLAN.EQ.l) ,0R.!NPLAN.EQ.4) ) GO TO 90 

THIS SECTION THEN TAKES THE SEWAGE INPUT THAT HAS BEEN GEN- 
ERATED AND PROCESSES IT ACCORDING TO ONE OF THE FIVE POLIC- 
IES. THE RESULTS FOR EACH HOUR ARE RECORDED IN A MATRIX. 

THE VALUE FOR THE MAXIMUM SEWAGE LEVEL IS ALSO RECORDED FOR 
EACH WEEKLY CYCLE 

ALPHA=ALPHA1 
90 DO 1000 J=l,24 ■ 

IF! ! J.EQ.l ) .AND. !I .EQ. 1) )G0 TO 855 
IF! ! J. EQ. 1) .AND. ! I .GE.2) )G0 TO 860 
Z!I,J)= Z! I,J-1)+X!I ,J)-ALPHA 
GO TO 880 

855 IF (NPLAN.EQ.4) ALPHA=ALPHA1 
Z( I, J)=X!I ,J)-ALPHA 
GO TO 880 

860 Z!I,J)= Z!I-1 ,24 ) +X! I , J ) -ALPHA 
880 IF!Z! I ,J) .LT.O.O)Z! I,J ) = 0.0 

IF!ZMAX!M,L) .LT.Z! I,J) ) ZMAX! M, L)=Z! I , J) 

TANK=Z!I,J) 

IF!NPLAN.EQ.l ) GO TO 1000 
IF!NPLAN.EQ.2) GO TO 1000 

IF !NPLAN.EQ.3)CALL PLAN3 ! TANK , AL PHA 1 , ALPHA , 8900 ) 

IF (NPLAN.EQ.4) CALL PLAN4!T ANK, AL PH A1 , ALPHA , ZSTOP , 
IZSTART ,8900) 

IF (NPLAN.EQ.5) CALL PLAN5! ALPHA 1 , ALPHA , KSTOP , 8900 ) 

900 CONTINUE 
1000 CONTINUE 
1500 CONTINUE 

THE PROCESSING TIME IS ALSO COMPUTED AND RECORDED 

IF (NPLAN.EQ.l) ITIME=168 
IF!NPLAN.EQ.2) ITIME=0 
IF (NPLAN.EQ.5) ITI ME= ( 24-KSTOP ) *7 
IPR0SS!M,L)=168-ITIME 



59 



SEWAGE LEVEL CAN BE PLOTTED FOR EACH 



* 



IF DESIRED THE WEEKLY 
HOUR 



* 

* 



IF(NODRAW.EQ.O) GO TO 2900 
N=0 

DO 2000 I=l»7 
DO 1900 J=l,24 
N=N+1 

YORD(N)=Z( It J) 

1900 CONTINUE 
2000 CONTINUE 

DO 2100 J=l»168 
2100 XORD(J)=J 

CALL DRAW! 168 , XORD , YORD , 0 t 0 » LABE L, TI TL E , 2 0 . 0 , 300 . 0 , 0 , 0 
ILAST ) 

* 3 ^ 

THE RESULTS CAN THEN ’BE PRINTED TO- THE EXTENT DESIRED 

2900 IF (NOPRIN.lt. 2) GO TO 2950 

WRITE{ 6»3500)ZMAX( M,L) ,ALPHA1,IPR0SS(M,L) . 

3500 FORMAT {' O' ,20X, • ZMAX=* tF8. 1, lOX, 'PROCESSING RATE = '' , 
1F6. 1 tlOX, ' PROCESSING TIME = '»I3t//) 

2950 IF(N0PRIN.LT.3) GO TO 5000 

WRITE(6»3000) { (Z( I » J ) » 1= 1 » 7 ) , J= 1 , 24 ) 

3000 FORMATCl Z MATRI X ' // { IX t 7 { F8 . 1 ,4X ) ) ) 

WRITE (6t 3100) { (X( I ,J ), 1=1,7) ,J=lt24) 

3100 FORMATCO X M ATR I X ' / / ( IX , 7 ( F8 . 2 , 3X ) ) ) 

5000 CONTINUE 
6000 CONTINUE 

THIS ROUTINE COMPUTES THE AVERAGE VALUES AND VARIANCES FOR 
THE MAXIMUM SEWAGE LEVEL AND THE WEEKLY PROCESSING TIME 

DIMENSION ZSUM(20) ,ZAVE(20) ,ZVSUM(20) ,TSUM(20) , 
1TAVE(20) ,TVSUM(20) ,ZVAR(20) ,TVAR(20) 

DO 7200 J=1,M 
ZSUM( J)=0. 

ZVSUM{ J )=0. 

TSUM( J)=0. 

TVSUM( J)=0. 

DO 7000 1=1, L 

ZSUM{ J)=ZSUM( J)+ZMAX( J, I ) 

TSUM{ J )=TSUM{ J)+IPROSS(J,I ) 

7000 CONTINUE 

ZAVE( J)=ZSUM(J)/L 
TAVE( J )=TSUM{J)/L 
DO 7100 1=1, L 

ZVSUM{ J)=ZVSUM(J)+{ (ZMAX(J, I)-ZAVE(J) )**2) 

TVSUM( J)=TVSUM( J )+{ ( IPROSSI J, I )-TA VE( J ) )**2) 

ZVAR{ J )=ZVSUM( J )/L 
TVAR( J ) =TVSUM( J) /L 
7100 CONTINUE 
7200 CONTINUE 

WRITE( 6,8000) NPLAN 

8000 FORMATCl ZMAX MATRIX FOR POLICY NO. ',12//) 

DO 8200 J=1,L 

WRITE! 6, 8100) { ZMAX { I , J) , 1=1 ,M) 

8100 FORMAT ( IX, 10(F8.2, 3X) ) 

8200 CONTINUE 

WRITE! 6,8500) !ZAVE! J) ,J=1,M) 

8500 FORMATCO AVERAGE ZMAX' / / ! 3X, 10 ! F 8. 2 , 2X ) ) ) 

WRITE (6,8600) ! ZVAR(J) , J=1,M) 

8600 FORMATCO ZMAX VARIANCE ' // ( 3X, 1 0 ! F 1 0 . 2 , 1 X ) ) ) 
WRITE!6,9000) 

9000 FORMATCl PROCESSING TIME MATRIX'//) 

DO 9200 J=1,L 

WRITE!6,8150 ) ! I PROSS ! I , J ) , I = 1 , M ) 

8150 FORMAT ( IX, 10! I 3, 3X) ) 

9200 CONTINUE 

WRITE (6,9300 ) (TAVE( J ) , J = 1,M ) 
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9300 FORMATCO AVERAGE T I ME ' / / ( 3X t 10 ( F6 . 2 , 2X ) ) ) 
WRITE (6, 9400 ) (TVAR(J ) , J=1,M) 

9400 FORMATCO TIME VARIANCE • // ( 3X, 10 ( F 10. 2 , IX) ) ) 
9999 CONTINUE 
STOP 
END 



FUNCTION GAMIALFA, BETA, START) 

C NON INTEGER GAMMA GENERATOR DEVELOPED BY D T PHILLIPS 
COMMON IX, I, J, ITIMEtNOPRIN 
IF( START.GT.l. 5) GO TO 58 
X3=1.0 

IF(ALFA.LE.2.0) GO TO 1 
IF( ALFA. LE. 5.0) GO TO 2 
GO TO 3 

1 B = 0.2177+(2.10299768*ALFA)-(4.34611961*ALFA’!==!=2) + 

. C(6.00914764’!'ALFA4'# 

C3)-(3.9509T728*ALFA=^=*4) + (0.97279251*ALFA=!=*5) 

GO TO 4 

2 B=0.64350+(0. 45839602* AL FA )-(0. 02952 801*ALFA**2 ) + 
C(0.00172718*ALFA* 

C*3)-(0. 00005810*ALFA**4)+( 0.00000082*ALFA**5) 

GO TO 4 

3 B=l. 33 408+ (0.22499991* ALFA )- ( 0 . 00230695*AL F A**2 ) + 
C(0.00001623*ALFA* 

C*3)-(0.00000006*ALFA**4) 

4 S=1.0+(1.0/B) 

START=5. 0 

12 IF(S-1 .0)110,60, 15 
15 S=S-1.0 
X3=X3*S 
GO TO 12 

110 GY=1. 0 + S* ( -0. 5771 0 17+S * (0. 985 854+S*(-0 . 87642 18+S*(0. 83 
1 S* (- 0. 5684729+ S* ( 0. 2 54 8205+ S*{ -0.0 51 49930) )))))) 
X3=X3*GY 

58 A=(X3/ (ALFA*BETA) )**B 
B=l. 0/B 
A=1.0/A 

60 CALI RANDUd X, IY,YFL) 

IX=IY 

GAM=(-A*ALOG(YFL) )**B 

RETURN 

END 



SUBROUTINE PLAN3 (T ANK , AL FAl , RAT E ,* ) 

THIS POLICY PROCESSES SEWAGE UNTIL 1600 AFTER WHICH TIME IT 
THEN EXAMINES THE LEVEL OF THE TANK AND IF BELOW 50 GAL. IT 
STOPS PROCESSING UNTIL 0500 THE FOLLOWING DAY. ON WEEKENDS 
PROCESSING STOPS ANYTIME THE LEVEL IN THE TANK REACHES 50 GA 

COMMON IX, I, J, ITIME,NOPRIN 
IF(( I.EQ.D.AND. (J.EQ.l) ) IFLAG = 0 
IF(N0PRIN.LT.2) GO TO 6 
IF((I.EQ.l).AND.(J.EQ.l) ) WRITE (6, 5) 

5 FORMATCO PROCESSING POLICY NO. 3 WAS USED ON THIS 
CRUN' ,//) 

6 IFCJ.EQ.D.AND.dFLAG.EQ.D) I T IM E= I T I ME + 24- JSTOP 
IF(J.EQ.l) IFLAG=0 

IF(I.LE.5) GO TO 10 
IF(TANK.LE.50. ) GO TO 30 

IFdRATE. LT. 0.01). AND. (TANK. GE. 500.) ) GO TO 20 
GO TO 40 

10 IF(J.LE.ll ) GO TO 40 

IF1TANK.LE.50. ) GO TO 30 
GO TO 40 
20 RATE=ALFA1 

IF (I FLAG.EQ.O) GO TO 40 
ITIME=ITIME+J-JSTOP 
IF(N0PRIN.LT.2) GO TO 26 
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WRITE( 6,25) I ,J 

25 FORMATC CONTINUED PROCESSING AT T= * , 1 1 , , 1 2 , // ) 

26 IFLAG=0 
GO TO 40 

30 RATE=0.0 

IF( IFLAG.EQ. DGO TO 40 
JSTOP=J 

IF(N0PRIN.LT.2) GO TO 36 
WRITE! 6,35) I ,J 

35 FORMATI'O STOPPED PROCESSING AT T= • , 1 1 , • - • , I 2 , / / ) 

36 IFLAG=1 

40 IF( ( I.EQ.7).AND.( J.EQ.24) ) GO TO 45 
GO TO 50 

45 IF(IFLAG.EQ.l) IT IME=ITIME+24-JST0P 
50 RETURN 1 
END 



SUBROUTINE PLAM4 (TANK, ALFA! , RATE , STOP, START, «) 

THIS IS A LEVEL CONTROL POLICY WHICH CONTINUOUSLY EXAMINES 
THE LEVEL OF SEWAGE IN THE TANK. IF BELOW SOME SET LEVEL 
(ZSTOP) THEN PROCESSING STOPS UNTIL THE SEWAGE LEVEL REACHES 
SOME UPPER VALUE (ZSTART) AT WHICH TIME PROCESSING CONTINUES 

COMMON IX, I, J, ITIME,MOPRIN 

IF{ ( I . EQ.l) .AND. ( J .EQ. 1) ) IFLAG=0 

IFINOPRIN.lt. 2) GO TO 6 

IF( ( I . EQ.l ) .AND. ( J.EQ. 1) ) WRITE! 6, 5) 

5 FORMATCO PROCESSING POLICY NO. 4 WAS USED ON THIS 
CRUN' ,//) 

6 IF! (J.EQ. 1) .AND.! IFLAG.EQ. 1) ) I T IM E= IT I ME + 24- JSTOP 
IF(J.EQ.l) JST0P=0 

IF(TANK.GT.STOP) GO TO 20 
RATE=0.0 

IF! IFLAG.EQ. DGO TO 40 
JSTOP=J 

IF!N0PRIN.LT.2) GO TO 11 
WRITE! 6,10) I , J 

10 FORMATCO STOPPED PROCESSING AT T= • , 1 1 , ' , 1 2 , // ) 

11 IFLAG=1 
GO TO 40 

20 IF!TANK.LE. START) GO TO 40 
RATE=ALFA1 

IF! IFLAG.EQ. 0) GO TO 40 
ITIME=ITIME+J-JSTOP 
IF!N0PRIN.LT.2) GO TO 35 
WRITE!6,30)I ,J 

30 FORMAT!' CONTINUED PROCESSING AT T= • , 1 1 , , I 2 , // ) 

35 IFLAG=0 

40 IF! ! I. EQ. 7). AND. (J.EQ.24) ) GO TO 45 
GO TO 50 

45 IF! IFLAG.EQ. 1) I T I ME= I T I ME+24- J STO P 
50 RETURNl 
END 



SUBROUTINE P LAN5 ! A LF A1 , RATE , KTI M E , * ) 

THIS IS A TIME CONTROL POLICY. IT ALLOWS PROCESSING FOR A 
FIXED PORTION OF THE DAY ONLY 

COMMON IX,I,J,ITIME,NOPRIN 

IF!N0PRIN.LT.2) GO TO 6 

IF! ! I . EQ.l) .AND. ! J .EQ. 1) ) WRITE! 6, 5) 

5 FORMATCO PROCESSING POLICY NO. 5 WAS USED ON THIS 
CRUN' ,//) 

6 IF! J.LE.KTIME) GO. TO 10 
RATE=0. 0 

10 RETURNl 
END 
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